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ABSTRACT 

The structure of the dusty circumstellar envelopes (CSEs) of proto-planetary nebulae (PPNs) reveals the 
mass-loss history of these sources and how such histories may differ for elliptical (SOLE) and bipolar (DU- 
PLEX) PPNs. To study the PPN structures via dust-scattered linearly polarized starlight, we have compiled the 
imaging -polarimetric data for all 18 evolved stars that have been obtained to date with NICMOS on-board the 
Hubble Space Telescope (HST). This alternative imaging technique provides a unique way to probe the dis- 
tribution of dust grains that scatter light around evolved stars. The new perspective gained from the imaging- 
polarimetric data has revealed several new aspects to the structures of PPNs. Point-symmetry is a prevalent 
imaging-polarimetric characteristic resulting from the azimuthal density gradient in the CSEs. Among these 
point-symmetric nebulae, three detailed morphological types can be differentiated by their polarized intensity, 
/ p , and polarization strength, P. While the azimuthal density gradient is reversed above and below the equatorial 
plane in optically thicker bipolar nebulae, there is no gradient reversal in optically thinner elliptical nebulae. 
The equatorial plane of the system defined by the integrated angle of polarization is not necessarily orthogonal 
to the axis of the apparent bipolar structure in the total intensity data. 

Subject headings: circumstellar matter — stars: AGB and post-AGB — stars: individual (IRAS z02229+6208, 
IRAS 04296+3429, IRAS 04395+3601, IRAS 06530-0213, IRAS 07134+1005, IRAS 
09452+1330, IRAS 10131+3049, IRAS 10178-5958, IRAS 10197-5750, IRAS 
12419-5414, IRAS 16594-4656, IRAS 17150-3224, IRAS 17245-3951, IRAS 
17423-1755, IRAS 17441-2411, IRAS 19343+2926, AFGL 2688) — stars: mass loss — 
reflection nebulae 



1. INTRODUCTION 

Proto-planetary nebulae (PPNs; e.g. iKwokl 11993b 
VanWinckel 2003) are low to intermediate initial mass 
(~ 1 — 8 Mq) objects that are rapidly becoming planetary 
nebulae (PNs) after evolving off the asymptotic giant branch 
(AGB). A PPN consists of a post-AGB central star and a 
circumstellar envelope (CSE) of gas and dust that is created 
by mass loss from the central star during the AGB phase. 
When mass loss ceases at the end of the AGB phase, the CSE 
physically detaches from the central star and keeps coasting 
away during the PPN phase. Thus, PPNs preserve the AGB 
mass-loss history imprinted in the density distribution of the 
CSEs. Therefore, PPNs are the best astrophysical laboratories 
in which we investigate the nature of dusty mass loss during 
the AGB phase. 

During the last decade, high-resolution imaging observa- 
tions have revealed a ple thora of morphologies among PNs 
and PPNs (see review bv lBalick & Frankll2002l) . One of the 
most important tasks in the AGB stellar research is therefore 
to determine the main physical factors that cause such a great 
morphological diversity emerging from AGB mass-loss ge- 
ometry, which is highly sp herically symmetric at the b egin- 
ning of the AGB phase (e.g. lHabing & Blommaertll993l) . The 
observed morphological diversity suggests a number of shell 
shaping processes, and various mechanisms have been pro- 
posed to explain the observed morphologies (Bal ick & Frank! 
2002, and references therein). In all possible cases, how- 
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ever, the CSE structure must deviate from the initial spherical 
symmetry. Therefore, there may exist a universal mechanism 
through which the required symmetry breaking is achieved, 
and thus, we have been performing surveys to identify the 
mechanisms that generate the initial non-sphericity in the 
CSEs of AGB stars. 

Using mid-IR (8-21/im) data. iMeixner et all (1 19991) have 
shown that all the extended PPNs in their sample are axisym- 
metric and fall into one of the two morphological classes: 
toroidal (showing limb-brightened peaks that suggest the 
presence of a dusty torus) or core/elliptical (showing an un- 
resolved emission core surrounded by a low surface bright- 
ness elliptical nebula). They concluded that the last stage 
of AGB mass loss is short-lived and inherently toroidal (i.e., 
equatorially-enhanced), probably coinciding with a period of 
enhanced mass loss known a s the "superwind" phase dRenzinil 
U98lHlben & Renzinilll983l) . 

In an optical imaging survey o f PPN candidate s using the 
Hubble Space Telescope (HST), lUeta et al.l d2000l) have de- 
tected elongated reflection nebulosities whose morphologies 
appear to bifurcate depending on the degree of the equatorial 
density enhancement in the CSE. SOLE (Star-Obvious Low- 
level Elongated) PPNs show a bright central star embedded 
in a faint, elliptically-elongated nebulosity, while DUPLEX 
(DUst-Prominent Longitudinally Extended) PPNs exhibit a 
remarkable bipolar structure with a dust lane completely or 
partially obscuring the central star. These optical morpholog- 
ical classes appear to possess a one-to-one correspondence 
to the mid-IR counterparts (core/elliptical w DUPLEX and 
toroidal w SOLE). 

Subsequent radia tive transfer mode li ng of some se- 
lected sources (e.g. [Meixner et al.1 119971: lUeta et al.l l2001bh 
IMeixner et~ai1 l2002, 2004 have shown that (1) the observed 
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PPN axisymmetry emerges during a brief stage of an en- 
hanced mass loss (~ 3 x 10~ 5 Mq yr" 1 ) towards the end of the 
AGB phase, (2) the enhanced mass-loss rate into the equato- 
rial plane is very much higher (;> 10 3 ) than that into the polar 
directions, and (3) the differences between SOLE and DU- 
PLEX PPNs stem from the physical distinction between the 
two types and cannot be attributed solely to the inclination an- 
gle effects. Furthermore, it appears that SOLE and DUPLEX 
PPNs develop respectively from low and high mass progen- 
itors. This agrees with a recent discovery of mostly bipolar 
structures among OH/IR-selected sources, which are located 
close to the galactic plane (i.e., more m assive), observed by 
HST dSahai & Sanchez Contrerasll2004l) . 

Imaging polarimetry is an alternative technique to probe 
the CSE structure, in which the dust-scattered linearly polar- 
ized component of light can be separated from t he unpolarized 
component of light. G ledhill and collaborators (Gledh iTT et alj 
l200ltlGledhiTll2005l) have found, through UKIRT imaging- 
polarimetric surveys of PPN candidates, that the axisymmet- 
ric structure of the PPN shells fall into one of the two classes, 
shell and bipolar, depending on the amount of dust in the CSE. 
These classes respectively correspond to SOLE and DUPLEX 
PPNs. The results of the previous mid-IR and optical surveys 
have thus been confirmed. 

We have performed near-IR imaging polarimetry with HST 
to investigate the dust distribution in SOLE PPNs by making 
use of the polarized intensity (7 p ) maps. Imaging polarime- 
try is especially powerful for this purpose, because the tech- 
nique isolates dust-scattered light arising from the CSE and 
permits us to probe the CSE structure even in the vicinity of 
the bright central star; CSE structure is often obscured in to- 
tal intensity (7) maps, which are affected by the enormous 
point-spread-function (PSF) of the central star. Since I p de- 
pends only on the Stokes Q and U, which are free from the 
unpolarized component by definition, it can be extracted even 
from the PSF-affected data. We have successfully demon- 
strated that (1) SOLE PPNs consist of a hollow spheroidal 
shell which has a built-in equatorial density enhancement and 
(2) the variable degree of the equatorial density enhancement 
(i.e., the optical depth of the CSE) would result in a variety 
of apparent PPN str ucture even amo ng this specific group of 
optically thin PPNs dUeta et al.ll2005l hereafter Paper I). 

In the present work, we extend our investigation of the CSE 
structure via I p into the realm of DUPLEX PPNs. We ap- 
ply the technique of imaging polarimetry in investigating the 
structure of the bipolar lobes, in which the optical depth is 
generally low enough for the method to work. To date, there 
have been observations of 18 evolved stars by HST/NICMOS 
in the polarimetric mode. Only a small portion of the entire 
data set has been published, and the I p data have rarely been 
utilized to study the CSE dust distribution. Therefore, we ana- 
lyze all the existing HST/NICMOS imaging polarimetric data, 
especially in I p , obtained from the archive in order to obtain 
more insights into the general mechanism(s) responsible for 
the structure development in the AGB wind shells. Below, 
we describe the data set and reduction procedure briefly (§2), 
present and discuss the results both individually and in groups 
(§3) and give a summary (§4). 

2. OBSERVATIONS AND DATA REDUCTION 

We have collected NICMOS imaging-polarimetric data for 
13 evolved stars from the HST archive and combined them 
with our own data for 4 SOLE PPNs that have already been 
published (Paper I). When PSF correction is necessary, we 



have used data for two PSF calibration sources, HD 12088 
and BD+32°3739, respectively in the NIC1 and NIC2 band. 
Table Q] summarizes the major observing parameters for all 
observations. Unfortunately, the data for IRAS 02168 and 
IRAS 10131 (in the POL-L band) are affected by severe sat- 
uration and not included in the subsequent analysis. In the 
following, we use the IRAS designation of the targets to refer 
to each object except for AFGL 2688, which does not have an 
IRAS designation. 

The standard set of NICMOS pipeline routines, CAL- 
NICA (Version 3.1), in IRAF/STSDAS 4 has been used in 
data reduction. Additional care has been taken to remove the 
pedestal effect (a variable quadrant bias) and the Mr. Stay- 
puft (the amplifier rin ging and streaking due to bright targets; 
iDickinson et al.ll2~002l) When mult iple frames are combin ed, 
we have used the drizzle package (Koe kemoer et alJl2002h in 
STSDAS. The resulting images maintain their original pixel 
scales of 0"043 pixel" 1 and 0."075 pixel" 1 respectively in the 
NIC 1 and NIC2 band. 

Derivation of th e Stokes parameters is done via the matri x 
inversion method dHines et alj|2000t IDickinson et al.ll2002l) . 
However, the coefficients are different for data taken before 
and after the servicing mission 3B (SM3B) in 2002 March, 
during which the NICMOS cooling system was installed. This 
resulted in two sets of the photometric conversion factors for 
the pre- and post-SM3B data. We have used these coefficients 
(listed in Table 2 of Paper I) where appropriate. A more de- 
tailed account of data reduction is found in Paper I. 

3. RESULTS AND DISCUSSION 
3.1. Results 

The imaging-polarimetric maps of evolved stars are pre- 
sented in a uniform format, following the scheme in Paper 
I (Figs. Q] to [14] Readers are encouraged to look at the color 
images available on the electronic version of the paper). In 
each figure, we show (a) the total intensity, /, map, (b) the po- 
larized intensity, I p , map, (c) the polarization strength, P, map 
and (d) the polarization vector position angle (PA.), 9, map, 
respectively from top to bottom. When the PSF-subtracted 
data are available, they are also displayed in the right column 
(panels e-h) in addition to the original data in the left pan- 
els, except for IRAS 10131 and AFGL 2688 for which the 
two-epoch data are presented (the earlier/later epoch on the 
left/right). 

The results are also summarized in Table [2] which lists the 
measured coordinates and / and I p fluxes of the object, the 
maximum P (P max ), the mean P {(P)) and its standard devi- 
ation, and description of the I p structure. The listed coordi- 
nates are the center of the polarization vector pattern that is 
determined by minimizing the sum of the square of the dis- 
tance from the presumed center to a line orthogonal to the 
polarization vectors. The pattern center is initially assumed 
to be the I peak (when seen), and the center finding process 
is iterated until the shift between the previous and current 
centers becomes smaller than the numerical accuracy of the 
analysis. The fluxes are determined by integrating the surface 
brightness over an aperture in which there is more than one 
fj s ky detection. P max and (P) plus its standard deviation are 
determined by using an aperture which is delineated by the 
well-established nebula edge or an aperture in which there is 

4 IRAF is distributed by the NOAO, which is operated by AURA under 
cooperative agreement with NSF. STSDAS is a product of STScI, which is 
operated by AURA for NASA 
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more than three er s k y detection. The I p structure is described 
by the overall shape, dimensions (typically major and minor 
axis lengths), and P.A. measured from the north to the east. 

To compare the present results with the previous unresolved 
polarization measurements, we have also performed aperture 
polarimetry and summarized the results with the measure- 
ments of Paper I in Table [3] We have done aperture photome- 
try on the three polarizer images and then used the measured 
fluxes from each polarizer image to compute the integrated 
P and 9. We have varied the aperture size from x 1 to x 8 
of FWHM in I to determine how the resulting values of P 
and 9 behave as a function of the aperture size. In general, P 
and 9 become stable and converge to a single value when the 
aperture size is equal to or larger than five times the FWHM 
size in /. Thus, the aperture size is defined to be the mini- 
mum aperture size (an integer multiple of FWHM in I) over 
which the integrated P and 9 are generally stable. The quoted 
errors in P (0.2 to 4.0%) and 6 (at most 6°) are thus the max- 
imum deviation in our measurements when the aperture size 
is varied beyond the adopted aperture size. These values are 
consistent with t he results of more thorough ana lyses done by 
iHines et al l J2000l) 5 and lBatcheldor et al.l d2Q06l) 6 . 

We have adopted the integrated 9 as a measure of the orien- 
tation of the equatorial plane because of the following reason. 
When the optical depth in the equatorial region is very high 
(r ^> 1), there is little I p because of the dilution effects due 
to multiple scattering. In this case, the integrated 9 tends to 
be parallel to the equatorial plane because of the orthogonal 
alignment of th e scattering plane with respect to the equato - 
rial plane (e.g., IWhitney & Hartmandfl99l lGledhillll2005h . 
However, when the optical depth is marginally high (~ 1), 6 
in the equatorial region can be orthogonal to the equatorial 
plane due to scattering at the walls of the bipolar lobes (e.g. , 
Johnson & Jonesl[T99ll ITrammell et al.lfT994t lGledhill200l 
Oppenheimer et al1l2005l) . Even in this case, the integrated 9 
tends to be parallel to the equatorial plane, because vectors 
in the lobes away from the equatorial region would still domi- 
nate. Therefore, the integrated 9 is more appropriate as a mea- 
sure of the orientation of the equatorial plane, and is expressed 
as # eq in the following. Note, however, that 9 eq is defined to 
be the angle perpendicular to the orientation of the equatorial 
plane so that the comparison between 6> eq and P.A. of the I p 
surface brightness distribution can be done intuitively. 

3.2. Individual Sources 
3.2.1. IRAS 04395+3601 (AFGL 618) 

IRAS 04395 (Fig. p is a DUPLEX PPN that is on the verge 
of becoming a PN (fWestbrook et al.lll975b ICalvet & Cohenl 
fl978h . The I map shows its well known bipolar structure 
elongated in the east-west dire ction due to multiple outflows 
(ITrammell. & Goodrichl 12002). The lobe interior structure 
is more clearly revealed in the I map t han in recent high - 
resolution ground-based near-IR maps (lUeta et al.l l2001al) . 
PSF subtraction did not improve the quality of the data be- 
cause no post-SM3B PSF data with a comparable degree of 
saturation is available and the whole PSF structure (especially 
of the spikes) is enlarged due to the bright, extended emission 
core. 

5 In verifying the matrix coefficients by measuring polarization of both po- 
larized and unpolarized standards, their measurements were consistent within 
0.3% in P and 5° in 

6 Possible instrumental polarization would make the matrix coefficients 
uncertain and was found to cause difference of 0.1% in P and 3° in 9. 



The 7 p map shows a general bipolar shape of the nebula 
superposed by clumpy structures in the lobes, a s has b een re- 
vealed in the optical by Tram mell. & Goodrichl (12002). How- 
ever, we do not see the tips of the outfl ows in I p . This i s 
consistent with their shock-excited nature (lUeta et al.ll2001al) . 
The morphology of the nebula changes dramatically in the P 
map. Strong polarization is localized in several clumps dis- 
tributed along the south edge of the west lobe and the north 
edge of the east lobe (respectively marked as W and E clumps 
in Fig. [lj). These clumps are distributed along a curve that 
delineates an inverse S. The 9 map shows general centrosym- 
metry, and the vector pattern center and the / peak are roughly 
coincident. 

The structure of the lobes in I p and P suggests that there 
is more dust on one side of the lobe than the other, given 
that the lobes are optically thin at near-IR. The regions 
of enhanced surface brightness in I p and P spatially corre- 
spond to the shortest of the multiple outflows seen in Ha 
(ITrammell. & Goodr ich 2002) 7 . This indicates that these re- 
gions of enhanced local density are associated with the lat- 
est outflows that channel recently ejected matter away from 
the central region, provided that mass loss proceeds at a con- 
stant velocity. In fact, 9 eq (83°) is aligned with the directions 
of these most recent outflows delineated by strong I p and P 
(~ 85°), and is not aligned with the bipolar axis (94°). This 
implies that the current orientation of the equatorial dust torus 
is oblique with respect to the bipolar axis (by about 10°). 

The dumpiness in the outflow structures of the enhanced I p 
and P surface brightnes s suggests mass loss that is episodic. 
ITrammell. & Goodrichl (|2002) have shown multiple ring-like 
structures in Ha and [S II] emission as well as an arc in Ha 
silhouette. These structures are in line with the episodic mass 
loss interpretation. These clumps, however, could alterna- 
tiv ely imply instabilities in the o utflows, as pointed out also 
bv lTrammell. & Goodrichl d2002h . 

Provided that outflows are associated with episodic mass 
ejection, the angular separation between clumps corresponds 
to the interval of mass ejection. The separations between 
clumps are 1"1 and l."3 in the west lobe respectively for 
the first and second closest clump pairs from the central star, 
and 0"9 in the east lobe. Using 120 km s -1 outflow veloc- 
ity (Trammell, & Goodrich 2002) and the distance of 1.7 kpc 
dBujarrabal et alJll988l) and assuming a linear outflow, these 
separations translate to 60 to 90 years. Given the uncertainty 
in the exact geometry of the outflows and of the clumps, these 
values are reasonably close to the val ue of 40 years (scaled to 
our ch oice of the distance) derived bv lTrammell. & Goodrichl 
(2002) based on the separation between the ring-like struc- 
tures seen in Ha and [S II]. 

3.2.2. IRAS 09452+1330 (IRC +10 216) 

IRAS 09452 (F ig. [H> is one of the brightest thermal IR 
sources in the sky dBecklin et al.ll 1 9690 . and hence, one of the 
best-studied C-rich AGB star. The CSE (of DUPLEX type) is 
seen as a fragmented bipolar nebula having fou r major emis- 
sion components in the optical and near-IR ( [Skinner et alJ 
1 1998 1 : lOsterbart et al.ll2000t iTuthill et al.l 120001: IWeigelt et alJ 
120021 iMurakawa et al.l l2005l) . 

The / map shows the main nebulosity that is stretched at 
a P.A. of 22°. We recognize three fragments in the north 

7 The east clumps correspond to the eastern outflow "d" (the northernmost 
of the three outflows), while the west clumps coincide with the shortest west- 
er n outflow (the southernmost o utflow in the west) as presented in Figure 1 
of Trammell, & Goodrich 1 2002). 
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(marked as the "NE", "N", and "NW" in Fig.|2^ respectively 
from east to west [PA. = 22°, -19°, & -66 1 ;Toste rbart et all 
120001: iTuthill et al J 12000: Wei gelt et al.ll2002h and an arc be- 
yond the south lobe (abo ut 3" away from the central peak; 
Mau ron & Hugginsll2000l) . 

The / p map shows the CSE structure more clearly than the 
I map, but it suffers from the polarization ghosts (high 7 p re- 
gions located radially with semi-reg ular intervals in th e west 
and southeast as indicated in Fig. [2j Hi nes et ail 120001) . The 
northeast and south lobes form the main bipolar nebulae at a 
P.A. of 22°. However, it is not clear if the north and north- 
west lobes also have their counterpart s in th e southeast be- 
cause of the ghosts. iMurakawa et al.1 (120051) have detected 
little (< 10%) polarization in the southeast in recent ground- 
based imaging polarimetry in the H band and suggested that 
the north and northwest lobes are caused by leakage of light 
through fissures in the west side of the central dust torus. The 
southeast ghosts are more obvious as in the P map. These 
ghosts are also the source of the anomalous polarization vec- 
tor pattern in the otherwise centrosymmetric 9 map. 

Our aperture polarimetry has yielded 9 eq of 32°. This is 
consistent with past polarimetric measurements in the opti- 
cal through mid-IR (IShawl & Zelmerfl 19701 iDvck et al.ll 197 ll 
Capps & Dvcklll972l:ICap ps & Knacke 1976: iTrammell etaU 
1994tlSerkowski & Shawll2001l) . Therefore, the bipolar lobes 
are obliquely oriented (10° off) with respect to the equato- 
rial plane. Moreover, the southwest lobe appears larger than 
the northeast lobe. This is probably due to inclination of the 
lobes with respect to the plane of the sky: the angle of 20° 
has been determined based on radiative transfer calculations 
dSkinner et al.lll998l) . 

3.2.3. IRAS 10131+3049 (CIT 6) 

IR AS 10131 (Fig. [3} is an IR-brig ht C-rich extrem e AGB 
star (lUlrich etalJll96fit ICoherifl979l LSchmidt et all (120021) 
have studied the CSE structure in detail using a set of HST 
data to show the presence of (1) large IR bipolar lobes (about 
10 times larger than in the optical) in the northeast-southwest 
direction and (2) multiple concentric arcs 1 " to 4" away from 
the central star. The two-epoch maps in Fig.[3]show the ghosts 
affecting all the quadrants of the field of view since the roll 
angle of HST at the two epochs differs by about 180°. We do 
not recognize any internal structure in the CSE core in the 7 P 
and P maps. 

A series of IR polarimetry between 1967 and 1976 
show e d that 9 varied between 110° and 180° dDvcket al.1 
I197U iKruszewskil |l97U IKruszewski & Covnd 119761). Es^ 
pecially, IKruszewski. & Covnd (119761) reported an IR polari- 
metric variabi l ity on the order of several months. Although 
ISchmidt et al.l ([2002) have determined that there is no evi- 
dence for temporal changes in the polarimetric properties of 
the CSE structure between the two epochs, there is some 6° 
change in 9 eq . This may be related to the intrinsic IR polari- 
metric variability of the object. With respect to the measured 
9 eq values, the shell elongation is not exactly perpendicular 
(8° to 14° offset). If 8 eq varies as 8, it may suggests variations 
in the dust density distribution near the equatorial plane on the 
order of months. 

3.2.4. IRAS 10178-5958 (Hen 3-401) 

IRAS 10178 (Fig.® is a PPN of DUPLEX type associated 
with a Be post-AGB star. HST observations in the optical 
and H2 1-0 S(l) bands have revealed this object's highly col - 
limated structure (iGarcfa-Lario et al.lll999tlSahai et all 1999). 



In the I map we see the collimated bipolar lobes that are 
linearly elongated into a P.A. of 73°. The equatorial region 
shows a pinched waist due to extinction. The inner structure 
in the IR does not appear as sharply as in the optical. 

The I p and P maps clearly reveal the physically thin nature 
of the walls of the cylindrical lobes. The dust lane that sepa- 
rates both lobes is more clearly delineated. Moreover, 7 p and 
P maps show, in the bipolar "skirt" region, that the northwest 
and southeast walls are more enhanced than the southwest and 
northwest walls, respectively. Furthermore, the southwest and 
northeast walls appear longer than the southeast and north- 
west walls, respectively. The 9 map indicates a rather aligned 
distribution of polarization angles by a shallow gradient of the 
image tone over the bipolar lobes (especially away from the 
central star). Nevertheless, the 9 map is generally centrosym- 
metric centered at the I peak. 

Our measurements yield 9 eq = 80°. The P map shows an el- 
liptically elongated region of weak P (<, 20% at a P.A. of 78°) 
superposed on the dust lane. This weak P region, therefore, is 
almost orthogonally oriented with the equatorial plane. This 
suggests that this weak P region may represent the region of 
diluted polarization due to starlight leaking through the bi- 
conical openings of the dust torus aligned with the equatorial 
plane. However, the cylindrical lobes are aligned at a P.A. 
of 73°, which is slanted by 7° with respect to the equatorial 
plane. In addition, the east end of the northwest wall and the 
west end of the southeast wall (the regions marked as "en- 
hanced skirt" in Fig. [4}:) seem to be slightly curved towards 
the central star with respect to the rest of the walls, which are 
very much linear. It appears as if these curved ends of the 
walls trace the weak P region. 

The polarimetric data thus uncover (1) the point-symmetric 
distribution of scattering matter in the cylindrical lobes (the 
northwest and southeast walls are more enhanced) and (2) the 
misalignment between the axes of the equatorial plane and 
the collimated cylindrical lobes. As a cause for the highly 
collimated structure of the nebula, Gar cfa-Lario et alJ (1999) 
on one hand suggested an extreme density contrast within 
the dust torus that would preferenti ally channel outflo ws into 
spatially-confined regions, while ISahai et all (1 19991) on the 
other proposed a dissipative hydrodynamic interaction of jets 
with the surroundings. Both scenarios, however, still re- 
quire additional mechanisms to direct the outflows/jets along 
the cylindrical lobes that are not orthogonally aligned with 
the equatorial plane. Alternatively, the enhancements in the 
northwest and southeast walls of the lobes may be simply due 
to the illumination effects as suggested from the shape of the 
weak P region, since the biconical openings of the central dust 
torus (implied by the orientation of the weak P region) are 
pointed towards these regions of enhanced I p and P. The pro- 
nounced double-shell structure at the east end of the northwest 
wall and the west end of the southeast wall may be related to 
the misalignment. 

3.2.5. IRAS 10197-5750 (Hen 3-404) 

IRAS 10197 (Fig. IS is a bipolar DUPLEX nebula sur- 
roundi ng an A2I star th at is associated with OH maser emis- 
sion (lAllen et alJ 1 1 980h . The optical reflection nebulosity 
has been imaged with HST to reveal its butterfly shape 
(ISahai etal J 11 9991). The HST imaging polarimetric data is 
generally consistent with the previous imaging polarimetric 
results presented by IScarrott & Scarrotl d 19951) . but reveal 
much more detailed structure. 

The / map exhibits some differences in comparison with the 
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optical map. Near-IR light delineates a figure S rather than a 
butterfly shape as in the optical maps, filling the regions of 
heavy optical extinction such as the equ atorial dust lane and 
the dark "spur' ' of dust dSahai et al.lll999h . The 7 p and P maps 
reveal the nebula's point-symmetry better than the 7 map: po- 
larized light is more concentrated on the east edge of the north 
lobe and the west edge of the south lobe. The distribution 
of polarized light emphasizes the figure S appearance of the 
nebula. The 9 map shows general centrosymmetry. However, 
the 7 peak position does not coincide with the center of the 
polarization pattern. The location of the central star derived 
from the polarization pattern is (a, (S)j2ooo = (10 h 21 m 33?88, 
-58°05'47."7). 

The 6>gq measurement yields 34°. Because of the complex- 
ity of the nebula, we can define a number of axes that may 
characterize the nebula structure as indicated in the 7 p map 
(Fig. The lobe axis is oriented with a P.A. of (1) 23° 
based on the overall structure of the polarized nebula, (2) -3° 
based on the tips of the figure S appearance, (3) 67° based on 
the strongest 7 p distribution in the nebula, and (4) -15° based 
on the protrusions along the western edge of the north lobe 
and the eastern edge of the south lobe. Furthermore, there is a 
region of weak P (< 20%) that is elliptically elongated at the 
center of the nebula (similar to what has been seen in IRAS 
10178). This weak P ellipse is oriented with a P.A. of 12° 
(indicated in the P map, Fig. [5};). 

None of these axes is exactly perpendicular to the equatorial 
plane. If we assume, as in IRAS 10178, that the orientation of 
the weak P ellipse corresponds to the direction of the biconi- 
cal openings of the central dust torus, then the torus is slanted 
with respect to the equatorial plane. In such a case, the cen- 
tral dust torus needs to precess with roughly 46° to shape the 
nebula structure by outflows through the biconical openings 
of the torus. 

3.2.6. IRAS 12419-5414 (Boomerang Nebula) 

IRAS 12419 is a very large optica l bipolar reflection neb- 
ula (55" x 21"; IWegner & Glasslll979l) . in which very strong 
optical polarization has been detected (45 - 60% in the lobes; 
iTavlor & Scarrottlll980l) . The HST data in Fig. |6]show polar- 
ization in the central 8" x 8" region of the nebula at 2{im. 
In the I map, we see an extremely faint nebulosity in the 
PSF-dominated central region, which appears to be shaped 
in a figure j at a P.A. of 0°. Although the nebula is too faint 
to perform PSF subtraction effectively, the / map shows that 
the near-IR lobes are elongated along the general north-south 
direction (slightly tilted along the northwest-southeast quad- 
rants). 

The 7 P and P maps marginally show the distribution of po- 
larized light along four "streaks" radially emanating from the 
central region (P.A. of -170°, -55°, 125°, and 10° as indi- 
cated in Fig. |6j>). These streaks seem to delineate the edges of 
the bipolar lobes. The 9 map pattern appears generally cen- 
trosymmetric. Although the degree of polarization is similar 
in the near-IR and optical around the central star (12 ±8% 
and 15%, respectively), the polarization angle is not the same 
(147° and 80°, respectively). This is probably caused by the 
matter in the equatorial region that is optically thick in the 
optical but optically thin in the near-IR. 

3.2.7. IRAS 16594-4656 

IRAS 16594 (Fig. |7]i is an evolved post-AGB star associ- 
ated with a multi-funneled optical reflect ion nebula, in which 
shocks start to shape the shell structure dHrivnak et al.lll999l; 



IVan de Steene. & van Hoof)l2003l) . The ima ging-polarim etric 
data of this object have been presented by lSu et all (120031) and 
us (Paper I). Since this object presents an interesting case of 
a rather optically thin DUPLEX PPN, we briefly re-describe 
the structure of the nebula. This object possesses (1) an el- 
liptical nebula elongated in the east-west direction (P.A. of 
82°) with an additional protrusion towards northwest (P.A. of 
-60°) and (2) the bipolar cusps in 7 p that delineate the in- 
nermost edges of the equatorial density enhancement in the 
shell. The cusps ar e spatially coincident with H2 emission 
dHrivnak et al.l20 04), which shows enhancements at the edges 
of the equatorial torus in an otherwise elliptically elongated 
shell (aligned with one of the optical funnels). P is slightly 
enhanced on the south edge in the east lobe and on the north 
edge in the west lobe. The vector pattern is generally cen- 
trosymmetric. 

9 eq is found to be 22°, which is unrelated to any of the 
observed structure. While IRAS 16594 is of DUPLEX type 
based on its shell structure, the fact that the central star is vis- 
ible in the optical and near-IR suggests that the optical depth 
in the equatorial region is near the lower limit for DUPLEX 
(about unity). Moreover, the two-pe aked mid-IR core struc- 
ture f ound in the equatorial region dGarcfa-Hernandez et al.1 
2004) also indicates rather low-density enhancement. Thus, 
our 9 eq measurement may have been affected by the bright 
PSF. 

3.2.8. IRAS 17150-3224 

IRAS 17150 (Fig.© is an archetyp ical DUPLEX PPN w ith 
a highly equatorially enhanced shell ( M eixner et al.ll2002l and 
references therein). The I map reveals the equatorial region 
of the nebula, which is very much obscured in the optical 
dKwok et al.lll998l: lUeta et alJl2000h. H ST images in near-IR 
bands have been presented bv lSu et all d2003l) . From the po- 
larization pattern, we find that the position of the central illu- 
mination source is (a, <5) J2 ooo = (10 h 17 m 19?81, -32°27'21."2), 
which is 0" 14 ± 0"07 off from the central I peak. 

The 7 p and P maps display the bipolar structure of the neb- 
ula, in which the hollowness of the lobes is highlighted by the 
polarized light delineating the edges and polar caps. In ad- 
dition, P is even more strengthened in the north wall of the 
northwest lobe and in the south wall of the southeast lobe. 
Thus, there is a preferential distribution of scattering matter 
in a figure S. The measured 6* eq is 120°, which is consistent 
with the bipolar axis at a P.A. of 122°. 

3.2.9. IRAS 17245-3951 

IRAS 17245 (Fig . @ is a less studied post-AGB star 
dHrivnak et alJ fl999l) . whose less evolved nature has been 
suggested by the star's association with OH maser emission 
dSevensterj|2002l) . The polarization maps reveal the bipolar 
nature of the nebula that does not posse ss an apparent bipo- 
lar nebulosity i n the total intensity maps ( Gledhil l et al.ll200"lT : 
ISu et alj|2003l) . While the main bipolar lobes are oriented 
along the north-northeast to south-southwest direction, the tail 
end of the nebula tips appears to be twisted towards oppo- 
site directions: the north tip to northwest and the south tip to 
southeast. This tendency is seen more clearly in the P map, 
in which the strongly polarized tips are found at a P.A. of 7°. 
The bipolar axis of the overall nebula elongation is oriented at 
a P.A. of 15°, while 9 eq is found to be 21°. Thus, there is 6° 
misalignment between the bipolar axis and the orientation of 
the equatorial plane. 
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3.2.10. IRAS 17423-1755 (Hen 3-1475) 

IRAS 17423 (Fig. [I0]l is a highly collimated bipolar PPN 
(of DUPLEX) associated with h i gh velocity outflow s (e.g. 
ISanchez Contreras & Sahal 1200 It iRiera et al.1 12003b - The 
present data cover only part of the brighter northern lobe in 
the NIC1 band. The / map mainly shows the northwest lobe of 
the central hourglass-shaped bipolar structure that is inclined 
with respect to the plane of the sky. The extended knots ar e 
seen in the / map (up to the NW2 knot; IRiera et al.ll2003l) . 
The 7 p map reveals the physically thin walls of the bipolar 
lobes more clearly. 

Strong polarization is found close to the west side of the 
north lobe and to the east side of the south lobe (in the central 
"hourglass" part). These strong P regions are found at a PA. 
of -54° (indicated in Fig.fTOb). This angle is nearly orthogo- 
nal to the equatorial plane (8 eq = 128°). There is 5° misalign- 
ment between the equatorial plane and the general orientation 
of the bipolar lobes (at a PA. of 133°). Thus, these strong 
P regions can be interpreted as the particular locations in the 
walls of the bipolar lobes that are illuminated by the central 
star through the biconical openings of the central dust torus, 
which is not orthogonal to the axis of the bipolar lobes. 

3.2.11. IRAS 17441-2411 

IRAS 17441 (Fig.[TT) is a G0I post-AGB star associated 
with an optical nebula of DUPLEX type (ISu et all 119981: 
lUeta et alll2000h . The figure S shape of the nebula is more 
pronounced in the near-IR than in the optical, as seen in the I 
map. The I peak coincides with the polarization pattern cen- 
ter, and thus, represents the central source obscured in the op- 
tical. 

As in IRAS 10197, a number of axes can be defined based 
on various morphological/polarimetric characteristics. The 
tips of the figure S shell point to a PA. of -2°, while the inner 
region of the nebula, where I p is the strongest, is elongated 
along a PA. of 46°. The P map reveals the walls of the bipo- 
lar lobes, which are generally oriented towards a PA. of 6°. 
However, the central weak P region (< 20%) appears elon- 
gated along a PA. of 20°. Aperture polarimetry determines 
eq of 28°, which agrees with none of the characteristic PA.s 
d escribed above. 

Oppe nheimer et all (120051) have constructed a detailed 
model of this object using evacuated bipolar lobes with an 
equatorial disk by fitting their polarization data from optical 
to near-IR. Their model is quite successful in reproducing the 
highly polarized nature of the shell as well as the axisymmet- 
ric overall geometry of the object. However, their model does 
not address the point-symmetric aspect of the shell structure. 
The polarization structure is peculiar in IRAS 17441 because 
I p is enhanced in the east edge of the north lobe and the west 
edge of the south lobe, but P is enhanced in the opposite side 
in each lobe. 

3.2.12. IRAS 19343+2926 (M 1-92) 

IRAS 19343 (Fig. Q2JI is an O-rich PN of DUPLEX type 
which has been extensively s tudied (iTrammell. & Goodrichl 
Il996t iBuiarraba Fet al.lll998alfR and references therein). The 
bipolar lobes that are separated by the dust lane are seen 
in the I maps as have been seen i n the previous qbserva- 
tions in atomic and H2 line em ission (Buian-ab aTeTai][T9 98b) 
as well as in CO dBujarrabal eTail 1l994l Il998al) . We do 
not see, however, the chains of spots in the bipolar cavities 
that have been identified as shock-excited clumps of mate- 
rial (|Bujarrabal et al.ll998bl) . Their non-detection in scattered 



light is consistent with their proposed shock-excited nature 
and suggests very little presence of scattering matter around 
the shocked gas. 

The overall polarization characteristics are generally 
consistent with the previous optical imaging polarimetry 
(ISchmidt et al.lll978l) . The present data show that (1) 7 p is 
concentrated in the walls of the bipolar lobes, revealing the 
hollow nature of the lobes, and (2) strong P is found along 
the north edge of the northwest lobe and the south edge of 
the southeast lobe of the nebula that is generally centrosym- 
metric. The nebula shape defines the symmetric axis (PA. of 
131°), along which CO outflow vel ocity is found to have a n 
almost constant latitudinal gradient (Buiarra bal et al.l [l998a). 
This axis is, however, not aligned with 9 eq of 15 1 °. As in other 
objects, we find some P enhancements in the nebula along the 
opposite edges in each lobe. 

3.2.13. AFGL 2688 (Egg Nebula ) 

AFGL 2688 (Figs. [O] and [Hi is a DUPLEX PPN that is 
consi dered one of the archetypes (iNev et al.l 1 9751; ISahai et all 
1 998), and its stru cture has recently been reviewed in detail by 
iGoto et all (120021) . The I maps show the object's well known 
bipolar lobes plus the concentric arc structure and the "spin- 
dle" structure along the dust lane between the lobes (seen only 
in the NIC2 band; Fig. fT4h). The I p maps display the point- 
symmetric nature of the lobes with the brightness enhance- 
ment appearing near the southeast corner of the north lobe 
and near the northwest corner of the south lobe. It is remark- 
able that the spindle structure has almost zero polarization. 
This is because the spindle structure is largely in H2 emission 
(ICox et al.lll997l:lKastner et al.ll200l . 

In the P maps, strong polarization arises near the north- 
west corner of the north lobe and near the southeast corner 
of the south lobe. The equatorial region of the nebula is op- 
tically thick even to near-IR light, but optically thin lobes 
show a generally centrosymmetric pattern in the 9 map from 
which the position of the illumination source is determined 
dWeintraub et al.ll2000l) . 6 eq (17° to 19° PA.) is slightly mis- 
aligned with respect t o the axis of the nebula (PA. of 12°; 
IWeintraub et al.ll2000t IGoto et al.f2002l) . Proper-motion mea- 
surements of the shell structure using the two-epoch NIC- 
MOS data have yielded the distance of 420 pc to the object 
(lUeta et al.ll2006h . 

3.2.14. IRAS 07134+1005, 06530-0213, 04296+3429, and 
(z)02229+6208 

Imaging-polarimetric characteristics of these SOLE objects 
have already been discussed in detail (Paper I). Here, we 
briefly summarize our findings. 

IR AS 07134 (Fig 3 of Paper I) is a C-rich F5Iab post-AGB 
star jHrivnak et al. 1989) associated with an elliptic al neb- 
ula dMeixner et alj|l997t iDaval et al.| [T99l Uura et al.| [2000a; 
lUeta et al.l 120001; iKwok et alj|2002l) . The shell is a slightly 
prolate spheroid with an inner cavity and a built-in equatorial 
enhancement, and hence, is not simply "toroid al" as has been 
interpreted by previous mid-IR observat i ons (iMeixner et al] 
U997l:lDayal et al.ll 19981: Uura et al.ll2000at IKwok et alJl2002h . 
In addition, the I p map confirmed that th ere is more matter o n 
the east side of the nebula than the west (Meixner et al.l2004T) . 

IR AS 06530 (Fig. 4 of Pap er I) is a C-rich F5Ib post-AGB 
star (Hrivnak & Reddy 2003) with an elliptical reflection neb- 
ula dUeta et al.l l2000). The nebula is a nearly edge-on prolate 
spheroid with a possible inner cavity. The nebula's density 
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structure appears such that the equatorial region of the cav- 
ity is filled with some matter, which manifests itself as the 
density-enhanced bipolar cusp-like structure. However, the 
total amount of dust along the equator is still not large enough 
to cause self-extinction within the shell (but close to the upper 
limit for a SOLE nebula). 

IR AS 04296 (Fig . 5 of Paper I) is a C-rich GOIa post-AGB 
star (lHrivnakl [l995) t hat is associated with an X-shape d opti- 
cal_elongafion dSahail [19991: lUeta et alj|2000b . While 
(1999) suggested that these elongations represent a pair of 
bipolar jets and a collimating disk, we interpreted that the 
shell had two, most likely hollow, spheroids: one (PA. of 26°) 
is nearly edge-on and the other (PA. of 99°) is inclined with 
respect to the plane of the sky. The two-spheroid interpreta- 
tion w as favored by a recent dus t scattering modeling of the 
object (Oprjenheimer et al. 2005). 

IRAS 02 229 (Fig. 6 of Paper I) is a C-rich G8-K0Ia post- 
AGB star (iHrivnak. & Kw ok 1999) su rrounded by an e llip- 
tical nebula of dust-scattered sta rlight (lUeta et al.ll2000i) and 
dust emission dKwok et al.ll2002l) . This object is the most pe- 
culiar among these four SOLE nebulae. This nebula does not 
show the hollow shell structure, but has a greater equatorial 
enhancement along the minor axis (PA. of 149°). These pe- 
culiarities have been attributed to a higher degree of equatorial 
enhancement in the shell and to the inclination angle of the 
nebula that is marginally edge-on. However, our 9 eq measure- 
ment (101°), which is not consistent with this interpretation, 
appears affected by the bright starlight in this marginally op- 
tically thin nebula. 

3.3. Discussion 
3.3.1. Imaging-Polarimetric Characteristics of D UP LEX PPNs 

Having reviewed the imaging-polarimetric characteristics 
of individual sources in §3.2, it now appears that the most 
prevalent characteristic among the observed DUPLEX PPNs 
is their point-symmetric appearance in the polarized surface 
brightness. This point-symmetric appearance of these nebu- 
lae in polarized light stems from the fact that one side of a 
lobe shows stronger polarization than the other and that the 
stronger polarization occurs at the opposite sides of the oppo- 
site lobes. 

Among 13 DUPLEX PPNs in this compilation, 10 sources 
exhibit a sign of point-symmetry while the degree of point- 
symmetry varies among the sources. The exceptions are IRAS 
09452, 10131, and 12419. In IRAS 10131 and 12419, the 
structure of the polarized surface brightness is not determined 
well since the polarized nebula is too compact in the former 
and the polarization flux is too weak in the latter. Thus, ne- 
glecting the two inconclusive cases, we have 10 out of 1 1 ob- 
jects that display point-symmetric imaging-polarimetric char- 
acteristics. In IRAS 09452, the main lobe pair (the NE and 
S lobes in Fig. |2]i is aligned along a straight line and so do 
the other lobes (the N and NW lobes) and there is no appar- 
ent indication of point-symmetric variation of polarization in 
the lobes. Thus, IRAS 09452 provides an example in which 
the imaging-polarimetric characteristics of the nebula are de- 
termined by the illumination effects. In fact, the N and NW 
lobes have been shown recently to be due to scattering of 
stray starlight leaking thro ugh fissures in the opt ically thick 
toroidal core of the CSE dMurakawa et all 12005). None of 
these 10 DUPLEX PPNs that possess point-symmetry, how- 
ever, shows obvious evidence for the illumination effects as 
does IRAS 09452. 



In DUPLEX PPNs, the equatorial regions can be very opti- 
cally thick in general as evidenced by the presence of a dust 
lane. Nevertheless, the bipolar lobes themselves are optically 
thin enough that they are seen via dust-scattered starlight. 
This allows us to adopt the same argument for SOLE PPNs 
that we used in Paper I and to investigate the structure of 
the lobes in DUPLEX PPNs by relating the 7 p and P surface 
brightnesses to the amount of dust along the line of sight. Fol- 
lowing a simple scattering angle effect (i.e., high degrees of 
polarization are expected from scattering at near 90° angles) 
as in Paper I, we can attribute the observed point-symmetric 
appearance in these 10 DUPLEX PPNs to the column density 
variations along the periphery of these nebulae. Since there 
is more dust on one side of a lobe than the other, we con- 
clude that there is an azimuthal density gradient in these lobes. 
Asymmetric illumination can potentially induce similar point- 
symmetric nebula morphologies in polarized flux. However, 
polarized flux enhancements can occur anywhere and do not 
have to occur at the periphery of the nebulae. This is incon- 
sistent with the observed trend in that the regions of enhanced 
polarization occurs near the periphery of the nebulae. Thus, 
we conclude that the asymmetric effect is less likely than the 
density effect as the cause for the point-symmetric appearance 
of the nebulae in polarized flux. 

While the overall imaging-polarimetric morphology of the 
majority of the sample sources can be described as point- 
symmetric, there appear to be three distinct CSE structures 
among those point-symmetric nebulae. There are nebulae that 
show polarized fluxes near the periphery of the lobes in both 
7 p and P (IRAS 10178, 16594, 17150, 17423, and 19343). 
As discussed for SOLE PPNs in Paper I, this structure sug- 
gests that the lobes are hollow shells. Such hollow lobes can 
be created by either the cessation of mass loss in the polar 
directions or the sweeping up of the mass-loss ejecta into 
the po lar directions. Recently, iGarcfa-Arrendondo & Frankl 
(2004) have shown hydro-dynamically that interactions be- 
tween an AGB wind and a jet from an accretion disk can cre- 
ate a narrow-waist bipolar nebula if the jet do minates. A fa- 
mous well-collimated bipolar PPN, M2-9 (e.g. Schwa rz"et al.l 
fl997h . which is very similar to IRAS 10178, has been inter- 
preted as a symbiotic system evolving into a PN based on a 
scenario in which ou tflows emanate from the accretion disk 
dLivio & Sokedl2001h . In fact, IRAS 19343 has been recently 
suggested to be a symbiotic system based on a detailed spec- 
troscopic study dArrieta et al.ll2.Q05b . Furthermore, the pres- 
ence of precessing jets has been suggested in IRAS 17423 
and 19343 based on th e presence of shock excited knots 
(iBuiarrabal et all Il998bt ISanchez Contreras & Sahall 120011: 
iRieraet al.ll200l . Thus, such systems may reasonably ex- 
plain the morphologies of these objects and even the oblique 
angle of the equatorial plane suspected in some of these ob- 
jects (see below) may well be explained under the presence of 
precessing jets. 

There are also nebulae that show P only near the periphery 
of the lobes but I p inside the lobes (IRAS 10197, 17245, and 
17441). These nebulae are also characterized by the very pro- 
nounced figure S shape. Unlike the previous case, the hollow 
P structure does not indicate the hollowness of the lobes since 
7 p is detected all over the lobes. The reduced P in the middle 
of the lobes is due to increased unpolarized light in the mid- 
dle of the lobes. Since the central star is not very visible in 
these sources, it is not likely that the central star contributes 
to the increased presence of unpolarized light. Thus, we in- 
terpret that the unpolarized component of light is increased 
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in the middle of these nebulae due to multiple scattering in 
their optically thicker equatorial regions. Such enhanced op- 
tical depths in the vicinity of the central star along the equa- 
torial plane can be attributed, for example, to the presence of 
a disk around the central star. It has been shown recently that 
a warped equatorial density enhancement can generate point- 
symmetric lobes by a two-fluid h ydro-dynamical simulations 
dIckell2003l:lRrikhorst etaT]|2004l) . Similarities between mod- 
els and the nebula structure (especially of IRAS 10197) are 
remarkable. 

There are yet another nebulae that display dumpiness in 
the polarized surface brightness distribution (IRAS 04395 and 
AFGL 2688). As previously noted, the clumpy structure of 
IRAS 04395 may be due to episodic mass loss or instabili- 
ties in the outflows. The bipolar lobes of AFGL 2688 also 
show signs of hollowness but their 7 p structure (especially of 
the northern lobe) is clumpy, and therefore, AFGL 2688 may 
belong to the previous group. However, none of the sources 
in the previous group shows the presence of the concentric 
arcs as in AFGL 2688. Thus, AFGL 2688 may be a very 
unique object. Although the nebula structure is inconclusive 
due to its compactness, IRAS 10131 may be related to AFGL 
2688 because IRAS 10131 also shows some arcs in the CSE 
dSchmidt et al.ll2002l) . 

3.3.2. SOLE and DUPLEX Morphologies 

In this section we compare the morphological characteris- 
tics of SOLE and DUPLEX PPNs in this the most complete 
compendium of imaging-polarimetric data for PPNs taken 
with HST including data presented in Paper I. In Paper I, we 
have shown that (1) optically thin SOLE PPNs consist of a 
hollow spheroidal shell with a built-in equatorial density en- 
hancement, which manifests itself as the toroidal structure and 
(2) there is a range of optical depths, which can influence the 
observed morphological traits even among SOLE CSEs. In 
the present work, we have concluded that there exist the az- 
imuthal density gradient in addition to equatorial density en- 
hancements in the bipolar lobes of optically thicker DUPLEX 
PPNs. 

The distinction between SOLE and DUPLEX PPNs was 
initially proposed based prima rily on the pres ence/absence of 
the optically thick dust lane dUeta et alj|2000l) . The present 
imaging polarimetry provides an alternative means to differ- 
entiate these two morphological types. As we have reviewed 
in $ 3.3.1, we have found a prevalent morphological feature 
- point-symmetry in the surface brightness distribution of po- 
larized flux. This point-symmetric appearance is probably due 
to the azimuthal density gradient being reversed in the oppos- 
ing lobes above and below the equatorial plane. On the other 
hand, in SOLE PPNs there is no reversal of the azimuthal den- 
sity gradient above and below the equatorial plane, although 
the azimuthal density gradient is present. The best example 
for this is IRAS 07 1 34 presented in Paper I, in which the pres- 
ence of the azimuthal density gradient has been shown in th e 
dust and CO gas distribution in its shell (jMeixner et al. 2004). 

Thus, it appears that mechanisms that generate a reversal 
in the azimuthal density gradient in the opposite lobes are 
related to the fundamental nature of the DUPLEX nebulae 
while these mechanisms do not operate in the SOLE nebulae. 
Therefore, we postulate that the emergence of the reversal in 
the azimuthal dust density distribution is directly linked to the 
origin of the morphological bifurcation between SOLE and 
DUPLEX PPNs. However, it is still inconclusive at this point 
what causes this reversal. 



3.3.3. Oblique Orientation of the Equatorial Plane 

In the present work, we have adopted 8 eq , the integrated 9 
over an aperture, as a measure of the orientation of the equa- 
torial plane (modulo 90° to be precise). In conventional imag- 
ing, the orientation of the equatorial plane is inferred from the 
appearance of the dust lane and/or of the bipolar lobes. As we 
have noted earlier, the way starlight gets scattered in the equa- 
torial regions can be strongly influenced by the local density 
enhancements. Therefore, the silhouette morphology of dust 
lanes, especially in high spatial resolution imaging data, can 
be susceptible to the local density distribution. In such cases, 
imaging polarimetry can alternatively constrain the orienta- 
tion of the equatorial plane better than the conventional imag- 
ing. 

We have seen in §3.3.2 that 6* eq does not necessarily corre- 
sponds to the orientation of the apparent dust lane seen in the I 
data. Among 13 DUPLEX PPNs, we have found in 10 sources 
that the discrepancy between 8 eq and the orientation of the 
bipolar axis is more than the quoted error, assuming that the 
uncertainty in the bipolar axis orientation is ~ 2°. The two an- 
gles are equal within errors in one case (IRAS 17150). In two 
cases 6> eq is not reliably fixed probably due to the bright central 
star and faint polarized surface brightnesses (IRAS 12419 and 
16594). In the 10 oblique cases, however, the determination 
of the orientation of the bipolar axis from the 7 p structure is 
somewhat subjective: the 2° uncertainty may be inappropriate 
for three cases (IRAS 17245, IRAS 17423, and AFGL 2688). 
Even so, there are still seven sources in which the equatorial 
plane is more than ~ 10° oblique from the perpendicular ori- 
entation with respect to the bipolar axis, and there are more 
than twice as many oblique cases as orthogonal cases. 

Even if we take a stricter stance on the subjectivity in defin- 
ing the bipolar axis and on our informal error analysis in de- 
riving 6> eq , there are cases in which the oblique angle is more 
than 20° (IRAS 17441 and 19343). This issue needs to be 
investigated further in future as it may be related to the mech- 
anisms to generate the equatorial density enhancement and 
bipolar lobes. 

4. SUMMARY 

We have presented a compendium of near-IR imaging- 
polarimetric data of evolved star CSEs that have been ob- 
tained with HST/NICMOS to date, in order to study polar- 
ization characteristics of these CSEs, especially of optically 
thicker bipolar DUPLEX PPNs. The data have shown that 
their point-symmetric appearance is prevalent (seen in 10 out 
of 1 1 objects). Since the prevalent point-symmetry does not 
seem to be due to the illumination effects in almost all cases 
(except for one), we have concluded that the observed point- 
symmetry is due to the presence of the azimuthal density gra- 
dient in the CSEs. 

These point-symmetric nebulae can be grouped into three 
based on the detailed morphologies. The first group con- 
sists of objects that show polarized fluxes near the periphery 
of the lobes in both 7 p and P (IRAS 10178, 16594, 17150, 
17423, and 19343). By using the same argument as in Paper 
I, we suggest that the bipolar lobes in these sources are hollow 
shells. The second group consists of objects that are charac- 
terized by the very pronounced figure S shape and show P 
only near the periphery of the lobes but I p inside the lobes 
(IRAS 10197, 17245, and 17441). Since the reduced P in the 
middle of the nebula is due to the increased unpolarized light 
in the nebula, we interpret that multiple scattering takes place 
in the relatively dense equatorial regions in these sources. The 
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last group is characterized by their clumpy nebulae (IRAS 
04395 and AFGL 2688) that may be due to episodic mass 
loss or instabilities in the outflows. While the lobes of AFGL 
2688 show similarities to the second group, the presence of 
the concentric arcs makes it rather unique. IRAS 10131 may 
be related to AFGL 2688 because of the presence of arcs. 

We have found that the ways in which the azimuthal den- 
sity gradient is manifested in SOLE and DUPLEX PPNs are 
very distinct: the direction of the azimuthal density gradient 
in DUPLEX PPNs is reversed in the opposing lobes above and 
below the equatorial plane whereas such a reversal is not seen 
in SOLE PPNs. The most probable explanation for this seems 
to be the presence of precessing outflows into the opposite 
directions and/or magnetic field that is reversed its polarity 
above and below the equatorial plane. We have also found 
indications that the equatorial plane of the system (defined by 
# e q) is not necessarily orthogonal to the axis of the apparent 
bipolar structure in the total intensity data (at least 7 out of 10 
sources). There are cases in which the oblique angle is more 
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TABLE 1 

HST/NICMOS Imaging -Pol arimetric Observations of Evolved Stars 















DITHER 




EXPTIME b 


ORIENTAL 






Source 


Alternative Names 


Date 


FILTER" 


SAMP_SEQ 


NSAMP 


PATTERN 


NPTS 


(sec) 


(deg) 


PID 


Reference 




Source Data 


IRAS 02168-0312 


Mira, O Cet 


1998 Aug 11 


POL-S 


STEP8 


13 


NONE 




135.8608 


34.7588 


7416 








1998 Aug 11 


POL-L 


STEP32 


4 


TWO-CHOP 


4 


169.8279 


33.9539 


7416 




IRAS (z)02229+6208 d 




2003 Mar 28 


POL-S 


STEP8 


16 


SPIRAL 


4 


319.6832 


169.035 


9377 


1 


IRAS 04296+3429 




2003 Mar 28 


POL-S 


STEP8 


19 


SPIRAL 


5 


519.5079 


-149.313 


9377 


1.2 


IRAS 04395+3601 


AFGL/CRL 618 


2003 Nov 26 


POL-L 


SPARS 64 


13 


SPIRAL 


4 


2559.868 


-35.4301 


9430 




IRAS 06530-0213 




2003 Mar 29 


POL-S 


STEP8 


19 


SPIRAL 


5 


519.5079 


-131.793 


9377 


1, 2 


IRAS 07134+1005 


HD 56126 


2003 Mar 29 


POL-S 


STEP8 


19 


SPIRAL 


5 


519.5079 


-129.432 


9377 


1 


IRAS 09452+1330 


IRC +10 216, CW Leo 


1997 Apr 30 


POL-S 


STEP8 


9 


NONE 




71.89689 


-113.462 


7120 


3,4 


IRAS 10131+3049 


CIT 6, RW LMi 


1998 May 5 


POL-S 


STEP 16 


12 


SPIRAL 


4 


319.7775 


-109.901 


7854 


5 
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STEP 16 


11 


SPIRAL 
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255.8052 


80.8889 
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5 
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MCAMRR 
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3.32561 
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3.02328 
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IRAS 10178-5958 
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124.507 


7840 
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Hen 3-404, Roberts 22 
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-133.434 
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-150.493 
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POL-L 


STEP 16 


24 


SPIRAL 
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815.5839 
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IRAS 17423-1755 


Hen 3-1475 


1997 Aug 9 


POL-S 
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21 
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239.7774 


-137.381 
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IRAS 17441-2411 




1998 Mar 10 


POL-L 


STEP32 


21 


SPIRAL 
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1247.867 


45.665 


7840 
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IRAS 19343+2926 


M 1-92 


1998 Apr 24 


POL-L 


STEP 128 


12 


XSTRIP 
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767.8638 


29.2312 


7839 




AFGL 2688 


Egg Nebula, V1610Cyg 


1998 Apr 8 


POL-S 


STEP 16 


20 


SPIRAL 


4 


831.556 


58.3121 


7423 


3 




2002 Oct 6 


POL-S 


STEP 16 


20 


SPIRAL 


4 


831.556 


-117.002 


9644 








1997 Apr 8 


POL-L 


STEP 16 


26 


SPIRAL 


4 


1215.39 


56.3451 


7115 


8, 9 






2002 Oct 6 


POL-L 


STEP 16 


26 


SPIRAL 


4 


1215.39 


-117.752 


9644 




PSF Calibration Source Data 


HD 12088 




2003 Mar 28 


POL-S 


STEP1 


19 


SPIRAL 


5 


79.7868 


162.51 


9377 


1 


BD +32° 3739 




1997 Sep 1 


POL-S 


SCAMRR 


9 


NONE 




11.368 


-92.7652 


7692 


1. 3 






2002 Sep 9 


POL-S 


SCAMRR 


9 


SPIRAL 


4 


12.992 


-108.88 


9644 


1 






1997 Sep 1 


POL-L 


STEP1 


9 


NONE 




41.86784 


-93.5733 


7692 


1. 3 






2003 Jun 8 


POL-L 


STEP1 


9 


SPIRAL 


4 


47.84896 


-109.63 


9644 


1 



REFERENCES. — l. lUeta et alj(2005h , 2. |GledhiH(20Ql , 3. |Hines et all (2003) , 4. |Osterbart eTai] (200Cl) , 5. lSchmidt et"al]f200l , 6. lSu et alj (20031) , 7. |Gledhill et"al](200l]) , 8. |Sahai et al]fl998l) , 9. |Weintraubetal](200Q|) 
d POL-S: short wavelength polarizers (0.8— 1.3/xrn, centered at 1. 1 /j,m); POL-L: long wavelength polarizers (1.89 — 2.1 /im, centered at 2.05 fim). Total exposure time per polarizer. c The ORIENTAT header parameter refers to PA of the image +y 
axis (degrees E of N).^ The "z" prefix in the IRAS designation is given to indicate the fact that this object was found in the Faint Source Reject File in the IRAS Faint Source Survey i Hrivnak et al. 1999). 



TABLE 2 

Summary of HST/NICMOS Imaging Polarimetry of the Circumstellar Shell around Evolved Stars 



Measured Coordinates (J2000)" Flux Polarized Surface Brightness Distribution 
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04 
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QA 
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3fSU 
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/ 


OS 
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JJ ± 14 


liLo + L.USp (nOllOWj 


9 7 v/ 1 n 
Z. / X 1 .u 


ZU 


IRAS 07134 


07 16 10.27 


+09 59 48.5 


0.08 ±0.27 


POL-S 


930() f 


2000 


81 


44 ±20 


ELS (hollow) 


4.8 x 4.0 


25 


TD A C (1(11 £ 9 

IKAS 0945 Z 


no /n ^^ 37 
09 4/ 5/.3/ 


ill 1 fi A 9 Q 
+ 13 10 42.5 
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a znn 
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65 
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Fragmented BPL 
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Ifl 1 fi n9 9fi 
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(1998-10-11) 


10 16 02.19 


+30 34 18.6 


0.08 ±0.16 


POL-S 


8470 


1540 


60 


25 =b 13 






35 


IRAS 10178 


10 19 32.48 


-60 13 28.8 


0.09 ± 0.07 


POL-L 


554 


146 


56 


27 ±12 


BPL (S) 


16.4x2.6 


73 


IRAS 10197 


10 21 33.88 


-58 05 47.7 


0.26 ±0.11 


POL-L 


1320 


330 


65 


29 ±14 


BPS (S) + Protrusions 


7.6x5.6 


23 


IRAS 12419 


12 44 46.07 


-54 31 13.5 


0.06 ±0.25 


POL-L 


380 


35 


41 


12± 8 


BPS (X) 


10.0 x 7.0 


163 


IRAS 16594 


17 03 10.04 


-47 00 27.0 


0.06 ±0.14 


POL-L 


620 


66 


73 


39 ±16 


ELS + Cusp 


5.0 X 2.2 


82 


IRAS 17150 


17 18 19.81 


-32 27 21.2 


0.14 ±0.07 


POL-L 


230 


60 


70 


26 ±18 


BPL (S) 


6.6x2.2 


122 


IRAS 17245 


17 28 04.58 


-39 53 45.1 


0.24 ±0.08 


POL-L 


145 


18 


48 


17±11 


BPL 


2.0 X 1.3 


15 


IRAS 17423 


17 45 14.23 


-17 56 47.0 


0.15 ±0.11 


POL-S 


2728 


488 


51 


26 ±12 


BPL (S) + Knots 


10.98 x 1.7 


133 


IRAS 17441 


17 47 13.51 


-24 12 52.0 


0.19±0.16 


POL-L 


212 


57 


59 


16±14 


BPL (S) 


4.0x2.1 


1 


IRAS 19343 


19 36 19.02 


+29 32 50.2 


0.27 ±0.22 


POL-L 


3360 


373 


52 


22 ±13 


BPL (S) 


12.9x3.9 


131 


AFGL 2688 (1998-04-08) 


21 02 18.66 


+36 41 37.4 


0.17 h 


POL-S 


4138 


2268 


89 


50±18 


BPL (inverse S) 


15.08 x 4.0 


12 


(2002-10-06) 


21 02 18.53 


+36 41 36.2 


0.23 h 


POL-S 


4698 


2728 


81 


51 zb 15 








(1997-04-08) 


21 02 18.76 


+36 41 36.9 


0.29 h 


POL-L 


5998 


2848 


74 


23 ±17 


BPL (inverse S) 


23.78 x 3.7 


12 


(2002-10-06) 


21 02 18.55 


+36 41 36.0 


0.22 h 


POL-L 


637« 


2958 


73 


23 ±16 









a Of the centrosymmetric polarization pattern center. Offset is the angular distance from the pattern center to the / peak with an associated error in the pattern center determination. POL-S: short wavelength polarizers (0.8 -\.3ftm, centered at 
l.ljum); POL-L: long wavelength polarizers (1.89-2.1^tm, centered at 2.05/im). c P max is the maximum polarization strength detected, while (P) refers the mean P and its standard deviation. 11 ELS: elliptical shell; BPS: bipolar shell without any 
apparent dust lane; BPL: bipolar lobes separated by a dust lane. e Of the nebula (counter-clockwise east of north).' The measured flux value is three times higher than the previously observed value (see §4 in Paper I). 8 The observed field did not 
cover the entire nebula. 11 No / peak. 
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TABLE 3 

Summary of Aperture Polarimetry 







Aperture 


pa 


9eq J 


P.A. b 


Source 


Band 


(arcsec) 


(%) 


(deg) 


(deg) 


TR AS 0999Q 


POL-S 


0.4 


6.1 ± 1.0 


101 ± 1 


59 


TR AS 049 Qfi 


POL-S 


0.6 


5.7 ±1.0 


1 1 ± 1 


26 99 




POL-L 


1 2 


/ . / m u. j 


Q'i -L- 1 


91 


TR AS OfWO 


POL-S 


1.6 


2.7 ±1.5 


159 ±2 


20 


TR A S 07 1 34 


POL-S 


6 


AO-Lin 


1 1£ 4- 1 

1 JO ZC 1 


L J 


TR AS 0Q4S9 


POL-S 


8 


1 1.2 ± 1.0 


32 ± 1 


22 


IRAS 10111 n998-05-051 


POL-S 


0.6 


4 6 _j_ 07 


49 ± 1 


35 


H 998-1 0-1 11 


POL-S 


0.6 


5.4 ±0.7 


43 ± 1 


35 


IRAS 10178 


POL-L 


1.0 


3.0±1.5 


80±1 


73 


IRAS 10197 


POL-L 


3.0 


8.5 ± 1.0 


34 c 


23 


IRAS 12419 


POL-L 


0.9 


3.5 ±1.0 


57 c 


163 


IRAS 16594 


POL-L 


0.9 


0.7 ±0.2 


22±6 


82 


IRAS 17150 


POL-L 


2.4 


17.5 ±2.0 


120 ±1 


122 


IRAS 17245 


POL-L 


1.1 


7.1 ±0.3 


21 c 


15 


IRAS 17423 


POL-S 


2.2 


5.2±0.1 


128±2 


133 


IRAS 17441 


POL-L 


1.0 


11.7 ± 1.0 


28±4 


1 


IRAS 19343 


POL-L 


0.9 


2.3 ±4.0 


151 c 


131 


AFGL 2688 (1997-04-08) 


POL-L 


3.8 


42.5 ±1.0 


19 z± 1 


12 


(2002-10-06) 


POL-L 


3.0 


37.3 ±0.8 


17 z± 1 


12 


a These errors are not via a formal error analysis. See te 


xt for details.' 3 Of the nebula (counter-clockwis 


eeast of north). Typieally ±2° en 


or. Reproduced from TableTj 1 / 3 Error less titan 0.5 deg. 
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5 0-5 
ARA (aresec) 

FIG. 1. — Imaging polarimetric maps of IRAS 04395+3601 in the (a) total intensity, /, (b) polarized intensity, 7 p , (c) polarization strength, P, and (d) 
polarization P.A., 9. The maps are in the standard orientation (north is up, east to the left) and centered at the 9 center with tickmarks showing the right 
ascension and declination offsets in arcseconds. The wedges indicate the scale of the image tone in data number per second (DN s , which is to be multiplied 
by PHOTOFNU [see Table 2 in Paper I] to convert to Jy) for / and / p , in percent for P, and degrees east of north in 9 (i.e. P.A. 0° means the polarization vector, 
which is perpendicular to the scattering plane, is oriented in the north-south direction). The east and west clumps are indicated. 




2 0-2 

ARA {arc sec) 

FIG. 2.— Same as Fig.[Q but for IRAS 09452+1330. Each lobe is identified with its designation, NE, N, NW, SE, and SW. 
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ARA (aresec) 

FIG. 3.— Same as Fig.Q] but for the two-epoch observations of IRAS 10131+3049 (Left: 1998 May 5, Right: 1998 Oct 11). The polarization ghosts are 
indicated. 
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4 2 -2 -4 

ARA (arc sec) 

FIG. 5. — Same as Fig.[T] but for IRAS 10197—5750. Radial line segments in the 7 p and P maps indicate the directions of possible axes of the shell. 




ARA (a re sec) 

FIG. 6. — Same as Fig.[T] but for IRAS 12419—5414. Radial line segments in the 7 p map indicate the directions of the bipolar walls. 
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2 -2 2 -3 

ARA (arcsec) 

FIG. 7.— Same as Fig.Q] but for IRAS 16594-4656 (Left: original, Right: PSF-subtracted). 




2 -2 
ARA (arcsec) 

FIG. 8.— Same as Fig.Q] but for IRAS 17150-3224. 
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1 0.5 -0.5-1 

ARA (arcsec) 

FIG. 9.— Same as Fig.Q but for IRAS 17245-3951. 
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2 0-2-4 

ARA (arcsec) 

FIG. 10. — Same as Fig.Q] but for IRAS 17423-1755. The position of the NWlaknot IRiera et all2003l) is shown in the / map. The radial line segment in the 
P map indicates the orientation of the strong P regions in the bipolar lobes. 
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Fig. 11.— 



ARA (arcsee) 

Same as Fig.[T] but for IRAS 17441— 241 1. Radial line segments in the P map indicate the directions of possible axes of the shell. 
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ARA (arcsec) 

FIG. 12.— Same as Fig.|7] but for IRAS 19343+2926 (Left: original, Right: PSF-subtracted). 
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5 -5 5 -5 
ARA (arcsec) 

FIG. 13.— Same as Fig.[j] but for the two-epoch observations of AFGL 2688 in the NIC1 band (POL-S; Left: 1998 Apr 8, Right: 2002 Oct 6). 




&RA (aresec) 

FIG. 14.— Same as Fig.[l3] but for the NIC2 band (POL-L) observations. 



